A vortex description of the first-order phase transition in type-I superconductors 
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Using both analytical arguments and detailed numerical evidence we show that the first order 
transition in the type-I 2D Abelian Higgs model can be understood in terms of the statistical 
mechanics of vortices, which behave in this regime as an ensemble of attractive particles. The well- 
known instabilities of such ensembles are shown to be connected to the process of phase nucleation. 
By characterizing the equation of state for the vortex ensemble we show that the temperature for the 
onset of a clustering instability is in qualitative agreement with the critical temperature. Below this 
point the vortex ensemble collapses to a single cluster, which is a non-extensive phase, and disappears 
in the absence of net topological charge. The vortex description provides a detailed mechanism for 
the first order transition, which applies at arbitrarily weak type-I and is gauge invariant, unlike the 
usual field-theoretic considerations, which rely on asymptotically large gauge coupling. 
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The role of topological excitations in the dynamics and 
thermodynamics of gauge field theories is a subject of 
great promise ranging from the understanding of vortex 
phases in superconductors, necessary for practical appli- 
cations, to the clarification of the mechanisms of charge 
confinement in non- Abelian gauge theories, such as quan- 
tum chromodynamics. 

Topological excitations are important as finite energy 
vehicles of disorder. Thus, phase transitions between a 
state of long range (e.g. magnetic) order and disorder can 
sometimes be understood as the proliferation of topologi- 
cal excitations, each bringing about disorder comparable 
to its size jjj . This is true in the XY model in two spatial 
dimensions (2D), which displays a Kosterlitz-Thouless 
(KT) transition to a disordered state due to vortex pair 
unbinding [g]. The second order transition in the 2D 
Ising model can also be formulated in terms of domain 
wall percolation. Furthermore there is evidence that the 
second order transition in the 3D XY universality class 
is associated with vortex string proliferation . 

In this letter we show how a first- order (discontinu- 
ous) phase transition in a simple gauge theory can be 
understood in terms of topological excitations and pro- 
vide detailed numerical evidence in support of this view. 
Our results suggest a dual vortex gas picture of the tran- 
sition, in analogy to the KT case. In the type-I Abelian 
gauge theory, however, vortices have attractive interac- 
tions, leading to characteristic metastability and collapse. 

For its simplicity and close relationship to the XY 
model we study the Abelian-Higgs (or Landau-Ginzburg) 
model in 2D. The Lagrangian density C is 
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dfi+ieA^. We focus on the regime where e 2 is larger than 
the scalar coupling A. This model describes the long- 
wavelength behavior of an ideal type-I superconductor. 

The standard argument [Q] for a first order transition 
in gauge+scalar field theories relies on large e, for which 
the gauge field is very massive and can be integrated out. 
This is only justifiable as n = vA/e — * 0, as it requires a 
separation of scales between 'heavy' gauge degrees of free- 
dom, which do not participate in the transition dynam- 
ics, and 'light' scalar field fluctuations. The result is a 
'free energy' F[<j>] that is both non-convex and, generally, 
gauge dependent. Nevertheless F[<j>] yields the correct 
qualitative picture for certain aspects of the transition. 

A description of the critical system in terms of gauge 
invariant degrees of freedom, for arbitrary k < 1, is there- 
fore desirable and may shed light on the mechanism of 
the transition. While vortices exist only as fluctuations at 
high temperature they become the only long lived mag- 
netic excitations of (|IJ) at low temperatures. Moreover, 
arbitrarily low energy excitations can be produced by the 
superposition of vortices and anti- vortices. These argu- 
ments suggest that vortices are relevant degrees of free- 
dom at criticality. In type-I, vortices attract each other 
independently of the sign of their quantized flux (topo- 
logical charge) || . Thus Abrikosov vortex lattices are not 
formed in applied magnetic fields. Instead type-I super- 
conductors form a non-extensive multi-winding vortex, 
restoring the normal phase at its core. Between the nor- 
mal and superconducting phases type-I superconductors 
exhibit a first-order transition. 

Vortices are radial, static classical solutions obeying 



where <p is a complex scalar field, F^ v — d^A^ — d v A^ is 
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in temporal gauge Aq = 0, where <f> = a(r)e m9 . 9 is the 
polar angle, n is an integer, Ag{r) = n/er — a(r)/r, with 
boundary conditions o~(r = 0) = a(r = 0) = 0, a(r — > 
oo) = u, a(r — > oo) — > n/e. For r larger than the core 
size the vortex behaves like a point source for massive 
scalar and magnetic fields. Then the vortex profiles can 
be written 

cr(r) = v- f(r); f(r) = a s v q K Q (m s r), (4) 

Ti 

Ag(r) = -a<3 v m Ki(mar), (5) 

er 

where mj = V~\v, m G = ev, as, ag are dimension- 
less constants, and Ki, modified Bessel functions of or- 
der i. The profiles correspond to Yukawa (massive) 
charges in 2D, in contrast to the familiar Coulomb loga- 
rithmic vortex solutions of the the 2D XY model. 

In the XY model the importance of topological charges 
to the phase transition is demonstrated by rewriting the 
partition function in terms of vortex degrees of freedom 
||. Unfortunately because the Abelian Higgs system is 
non-Gaussian it is impossible to perform an exact dual 
transformation to a partition function written exclusively 
in terms of a one and two body vortex terms. It is never- 
theless possible to perform this re-writing approximately. 

We begin with a superposition ansatz for an arbitrary 
number N of vortices, by constructing scalar and gauge 
vortex fields centered at N different loci r^, i 6 {1, N} 
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Ai(r,n,...r N ) = Ai(\r - n\) + . . . + A N (\r - r N \). (7) 

This ansatz is exact when the vortices are all coincident 
or all widely separated. Substituting (^-|^) into the static 
part of the Hamiltonian gives 



H = ^2e l + [mi(x)V G (\x - y\)m j (y) 



(8) 
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where the terms not shown correspond to 3 and 4- 
body effects, which are negligible in a low-density vor- 
tex ensemble. The charges mi(x) — ±nS(x) are in- 
tegers of either sign, corresponding to quanta of mag- 
netic flux, whereas qt — \rrn\ is always positive 
For well-separated vortices the potentials are Vc{r) = 
a G v% q K (m G r), V s (r) = -a s v^ q K (m s r) @. v cq is 
a measure of cr(T) and as, a G are weakly varying with 
the couplings and have been computed numerically by 
Speight ||: for k < 1, a s > a G . Then the 2-body poten- 
tial for a pair of like-charge vortices is 



V(r) ~ -qiqjV 2 eq [a 2 s K (m s r) - a G K Q (m G r)] 



(9) 



demonstrating that in type-I, when k < 1, the scalar (at- 
tractive) part of the potential dominates the interaction. 
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FIG. 1. a) A characteristic time evolution of at 
T — 0.01844, close to phase coexistence, a jumps abruptly 
between the normal and superconducting phases, b) The to- 
tal vortex density p for the same trajectory as in a), p is a 
disorder parameter vanishing in the superconducting phase. 

The statistical mechanics of particles with a finite- 
range, soft-core attractive potential was studied in sev- 
eral simple settings ||, aimed at elucidating thermo- 
dynamic gravitational instabilities. Unfortunately the 
Yukawa gas of (0), was not among these. Nevertheless 
this class of systems share important qualitative proper- 
ties; namely they always display a first-order transition 
between an (almost) ideal gas state at high T and a non- 
extensive clustered phase at low T ||. The latter is not 
a thermodynamic phase. It consists of a single bound 
cluster containing most of the particles. Qualitatively 
this clustering transition occurs at T c i, such that the ki- 
netic energy equals the interaction energy per particle. 
This gives a rough estimate of the vortex gas clustering 
temperature H, T c \ ~ e/2 = 0.019 where e measures 
the strength of the interaction. We estimate e from [jlO| , 
by averaging the strength of vortex-vortex and vortex- 
antivortex pair interactions, e = (e vv + e V y)/2, where 
e vv ~ 0.62wg q and e„ e = 3.80Wp q . Vortex interactions are 
softened by a small order parameter, v cq = c(T c + ) — 0.13. 

The equation of state for the almost ideal (low density) 
gas of vortices at high T can be computed by standard 
cluster expansion methods. For simplicity we model the 
attractive potential as a square well with a strength e 
and an interaction length I. Then 
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where we neglected terms proportional to p n , n > 3. The 
correction to the ideal gas behavior is negative as ex- 
pected for an attractive potential. For high T, B2 van- 
ishes. Interactions are most important at low T and the 
pressure P vanishes at T c \: 

T ci ~ , 6 ; , . (11) 



In 



Using pnl 2 ~ 0.37 = &{T~) — er(T+), and e as above we 
obtain T c \ ~ 0.020. Both estimates of T c \ are compatible 
with the measured T c and coincide in the limit e/T <C 1 
and pnl 2 ~ 1. 




FIG. 2. a) The probability distribution of a at 

T — 0.01844, in the critical region, showing coexistence of 
the normal and superconducting phases, b) A hysteresis loop 
for the vortex density p obtained by cooling (triangles) and 
heating (circles) the system through the critical region. 

These estimates of T c \ are very qualitative in nature. 
For our choice of parameters the exact 2-vortex potential 
is not known. Moreover the virial expansion is notori- 
ously unreliable near a transition, especially in the pres- 
ence of clustering. We plan to examine the quantitative 
behavior of the vortex particle ensemble, through a direct 
thermodynamic study using either the exact intervortex 
potential or the Yukawa gas interaction of (ph. 

Nevertheless, a vortex ensemble picture leads to several 
qualitative predictions i) p will behave as a disorder pa- 
rameter, vanishing discontinuously in the superconduct- 
ing phase, ii) the vortex ensemble must show signs of 
a clustering instability in the metastable phase and iii) 
multi-charged vortices must be visible under sudden non- 
equilibrium cooling. 



To test these predictions we consider the field evo- 
lution in contact with a heat bath, given by a system 
of Langevin field equations. Gauge invariance demands 
that the evolution preserves Gauss' law. This constraint 
still allows for several classes of dynamical equations |l2) , 
characterized by different gauge invariant stochastic gen- 
erators. We choose the simple set {|<^>| 2 ,-E}, leading to 
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9t4>i = 7T 4 , 

d t Ei = (V x B)i + J i} J t = -e 2 \<f>\ 2 Ai - ee ab (p a di<f> b , 
d t A z = Ei+TjgdtEi+Tg, 

with Ei — dtAi , B = V x A. The details of this choice are 
irrelevant to the state of canonical thermal equilibrium 
reached at long times. The indices a, b refer to the 2 real 
components of <j), whereas i is a spatial vector index, is 
the totally anti-symmetric rank 2 tensor. The stochastic 
sources T obey fluctuation-dissipation relations 

(r Sig (x,t)r Stg (x',t , )) = 2r) Sig T5(x-x')5(t-t'), (13) 



with (r. 



0. We choose e = 1.5, A = 0.1, and solve 



a lattice (non-compact) version of Eq. (|l2|), with r) s 
ri g = 0.05, dt = 0.02 and dx = 0.5. 
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FIG. 3. The density of like-sign vortices p(r) in a disc of 
radius r around a vortex, normalized by the average vortex 
density. Plots correspond to the metastable phase (solid), 
the near-ideal gas phase at T > T c (short dashed) and to a 
random distribution (long dashed), with the same density as 
the metastable phase. The latter exhibits correlations typical 
of the subcritical percolation problem. Vortex clustering is 
maximal in the metastable phase, a premonitory sign of the 
dynamical clustering instability, and disappears for T>T C . 



Fig. |l] shows that a substantial p ^ exists in the 
symmetric (normal) phase, but that all vortices suddenly 
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disappear when the system transits to the superconduct- 
ing state, as shown by the spatial average of a = \<fi\. 
Vortices in the normal phase are identified by their quan- 
tized fluxes, a quantity that is manifestly gauge invariant. 
The probability distribution of a close to T c is shown in 
Fig. [2^,. The double peak demonstrates coexistence of 
the two phases, characteristic of first order transitions. 
Fig. ||b shows a hysteresis loop in p, obtained by slowly 
heating and then cooling through the critical region. 




FIG. 4. Contours of magnetic flux after a fast temperature 
quench, showing the clustering of singly quantized vortices 
(smallest circular features) into large integer charge bound 
states. White (black), localized, regions denote (anti) vortices. 
The total collapse of the vortex ensemble was avoided due to 
fast cooling, which evaded the metastable region. 



The metastability of the vortex ensemble is a conse- 
quence of the small probability for the formation of a 
large vortex cluster. While the free energy is lowered 
through attractive interactions in the volume, the spatial 
cluster boundary, where vortices are rarefied, is thermo- 
dynamically costly. Thus small clusters are subcritical 
and can exist in the metastable phase without leading to 
its collapse. Evidence for incipient clustering is shown 
in Fig. where we plot the radial density of like-charge 
vortices around another vortex. Clustering is maximal in 
the metastable region, and negligible at higher T>T C . 

While we argue that vortices are the relevant degrees 
of freedom at criticality, their profiles cannot be clearly 
observed in the normal phase, because vortices appear 
there only as transient fluctuations. Below T c , where 
they could exist as well defined objects, their ensemble 
collapses and vortices disappear in the absence of quan- 
tized net flux. Fast quenches evading equilibrium in the 
metastable region do display well defined vortices and 
show striking evidence of their clustering, see Fig. ^. 

In conclusion, we argued for a vortex description of the 
mechanism underlying the first- order transition in the 2D 
type-I Abelian Higgs model and provided detailed numer- 
ical evidence in support of this picture. Below T c \ the 
vortex ensemble becomes metastable and eventually col- 



lapses to a non-extensive thermodynamic phase and the 
system becomes a superconductor. The vortex interpre- 
tation of the transition is gauge invariant and does not 
require e 2 ^> A, unlike the field-theoretic arguments for 
a first-order transition. Instead, the attractive nature of 
the vortex potential is manifest even in the weakest type- 
I regime and metastability and collapse are inescapable. 

The vortex description must now be subjected to closer 
quantitative scrutiny by direct studies of particle ensem- 
bles. Still, we believe that its qualitative success bodes 



well for applications to 3D |13|, where point vortices be- 
come lines, and may participate in interesting critical 
phenomena such as crystal melting and in cosmology. 
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